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Abstract 

Schrodinger 's equation with scalar and vector potentials is shown 
to describe "nothing but" hopping of a quantum particle on a lattice; 
any spatial variation of the hopping amplitudes acts like an external 
electric and/or magnetic field. The main point of the argument is 
the superposition principle for state vectors; Lagrangians, path inte- 
grals, or classical Hamiltonians are not (!) required. Analogously, the 
Hamiltonian of the free electromagnetic field is obtained as a twofold 
continuum limit of unitary hopping in Z(N) link configuration space, 
if gauge invar iance and C and V symmetries are imposed. 

PACS: 03.65.Bz,03.70.+k,11.15.Ha 



1 Introduction 



For at least two standard quantum systems, canonical quantization (or other 
classical-to-quantum substitution rules) can be avoided; it can be replaced 
by an intrinsic quantum mechanical consideration of "hopping" in a discrete 
configuration space. This only requires to interpret a familiar tool from model 
building — hopping amplitudes — as a first-principle concept. 

Hopping amplitudes have a long tradition particularly in solid-state the- 
ory ^ . On a fundamental level they have been used in lattice gauge theory 
0, ^ for discretizing (not avoiding) path-integral actions. More recently, in 
the field of quantum computation, hopping parameters are being used as 
collision constants in unitary cellular automata ^ designed for efficient 
simulation of the Schrodinger equation ||^ or 1-photon and Weyl equation 
0. These latter applications differ in a crucial way from the viewpoint taken 
here, by assuming locality in conjunction with a finite, irreducible time step. 
It has proven to be a major challenge to design algorithms satisfying that 
computational requirement. Apart from technical complications, however, 
unitary cellular automata in some cases require configuration spaces larger 
than the physical ones. For example, local hopping rules in d spatial dimen- 
sions are found to require 2(i-component wave functions 0. Consequently, 
a real spinless particle (as opposed to its computer simulation) can have a 
unitary and local equation of motion only with respect to continuous time. 

Hopping amplitudes can do more than approximate or discretize processes 
originally defined otherwise. They necessarily emerge as coefficients of a 
superposition when a particle is prepared in a position eigenstate. The crucial 
axiom here is that the state of a quantum particle is completely specified by 
a position at one instant of time. To illustrate the idea, consider a particle 
confined to a 1-dimensional array of discrete positions at a spacing a. Let us 
work in the Heisenberg picture and denote by |n, t) the eigenstate of position 
X = na a.t time t. 

To prepare a position n at time t means to prepare a state with an 
uncertain position at time t + dt, because any motional information is lacking 
from \n,t). For dt small enough, the uncertainty only relates to positions n, 
n + 1, and n — 1. Furthermore, n + 1 and n — 1 will occur symmetrically if 
we assume the symmetries of a free particle. Thus 

|n, t) = a\n, t + dt) + (3\n + 1, t + dt) + (3\n - 1, t + dt) (1) 

where a and (3 are some numbers dependent on the size of the time step. For 
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dt — > we must have a ^ 1 and (3^0, hence 

a = 1 + aidt + C(dt2) (3 = (3idt + 0{dt^) 

Thus the basic hopping equation (|1|) converges to the differential equation 

= ai\n,t) + I3i\n + l,t) + I3i\n - l,t) (2) 
We now use the statistical interpretation of the scalar product. From 

{n,t\n',t) = 6n,n' 

we find by differentiating with respect to t and using that the coefficients 
ai and /5i must be purely imaginary. Finally, we consider a general state 
vector in the Heisenberg picture, 

\i;) = Y.^P{n,t)\n,t) (3) 

n 

We take d/dt, use (0), put x = na, and reexpress ai and /?i by 



U = {ai + 2Pi)ih — 



1 a^p 



2m ih 
Thus we find 

Ti^ il){x + a, t) + il){x — a,t) — 2tjj{x, t) 



ih—ip{x,t) = U^{x,t) 

dt 2m a 
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This equation converges to the free Schrodinger equation in the continuum 
limit a — > 0. 

In Section 2, the hopping-parameter description of a Schrodinger particle 
is discussed in full generality. Hopping amplitudes will not be restricted to 
next neighbours, and it will only be assumed that the hopping amplitudes 
realise the full translational and cubical symmetries of the lattice in (9(l/a^) 
while any inhomogeneities in the hopping process are at most of 0{l/a). 
Then a (trivial) renormalization scheme exists for the continuum limit a 
which leads to the standard nonrelativistic Schrodinger equation, with a 
vector potential and a scalar potential. 

In Section 3, the hopping-parameter approach is applied to quantum elec- 
trodynamics without charges and currents. This requires the discretization 
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of both, the values of a field u{x) and its spatial variable x. The reader of 
section 3 is assumed to be somewhat familiar with lattice gauge theory 
In fact, the model considered in this section is a Hamiltonian version of the 
intensively studied Z{N) lattice gauge theory The Hamilton operator 
of the electromagnetic field is recovered in the twofold limit of iV ^ oo and 
zero lattice spacing. Section 4 contains some concluding remarks. 



2 Schrodinger particle in 3 dimensions 

Consider a simple cubic lattice where x = an is the position vector of a site, 
a is the lattice spacing, and n an integer vector. The most general hopping 
equation for a single-component wave function as defined in (H) is 

ih—^lj{x,t) = '^K{x,n,t)'ijj{x + an,t) (4) 

n 

The factor of ih is only cosmetic, since the hopping parameters k{x, n) can 
be any complex numbers, so far. Conservation of probability requires 

k{x — an,n,t) = k{x, —n,t) (5) 

An important case of reference is that of a free particle, characterized by 
hopping parameters with the full symmetry of the lattice. Then k{x, n, t) = 
Ko(n) because of translational invariances. Cubic symmetry implies 

Ko{n) = Ko{-n) (6) 

so that all ko(^) are real because of (^). Most importantly, the symmetry 
also implies J2n i^o (n) riirij oc 5ij . A convenient parametrization is 

V Ko(n) niUj = 6ij (7) 

n 

The reduced parameter m will be identified as the particle mass later on; the 
sign of m is discussed in the Conclusions. In general, the sum in equation 
d^) need not converge. Assuming convergence here is the basis for the non- 
relativistic physics as it emerges in the form of the Schrodinger equation in 
the continuum limit. 
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To recover the Schrodinger equation, we Taylor-expand the displaced 
wave functions on the RHS of (^), 

ip{x + an, t) = ipi^j t) + aniViipix, t) + ^a^ni'rij'ViVjipix, t) + 0{a^) (8) 

Again, let us consider a free particle first. Inserting K,{x,n,t) = K,o{n) in 
and using (H), (^, and (0) we find 

d - - 

ih— ip{x, t) = Eoi^ix, t) - —V ■ VV^(f, t) + 0{a) (9) 

where Eq = J2n ^o(^) is certainly infinite but does not affect the shape of the 
wavefunctions. In contrast, the parameter m determines the particle mass 
and must be finite, as anticipated in definition 

Now we "turn on" deviations of the hopping parameters from no{n). Let 
us put 

ti{x,n,t) = KQ^n) + Ki{x,n,t) (10) 

Again, we insert (H) in (^). The multiplicative terms on the RHS of 
now are EQip{x,t) +J2n'^i{^^^yt) i^i^yt). The inhomogeneous term can be 
rewritten as ^ 

+ -n,t)) ip{x,t) 

n 

Using (^) and expanding the ensuing displaced argument, we obtain the 
following form of the multiplication operator: 

- ^ (^Ki{x, n, t) — Ki{x, n, t)^ + 2 X! '^^ ■ ft, t) + 0{a^Ki) 

n n 

This shows that for a finite if-dependent contribution, the real part of ki 
must be of 0{1) while the imaginary part can be of 0{l/a). Hence, if we 
define a vector potential 

A{x, t) = -y- n {x, n, t) (11) 

n 

then the multiplicative terms of (^ take the form 

(^Eo + J2 ^^li^, n, t)^ ^{x, t) +t^(v- A{x, t)) ^{x, t) (12) 



5 



The gradient terms on the RHS of (^) can be written as 



-V-ip{x, t) ■'^n {k{x, n, t) — k{x, —n, t)) 



By (p!OD and (13) this is equal to 



t) ■'^n (ki{x, n, t) — Ki(x — an, n, t 



The displacement of x in ki{x — an, n, t) produces a term of higher order in 
a which can be neglected in the limit a — 0. Thus the only relevant con- 
tribution to the gradient terms comes from the imaginary part of ki {x, n, t) 
and is of the form 

i^(yilj{x,t))- A{x,t) (13) 

where A{x, t) is the same as in ([TT|). 

With inhomogeneities of 0{1) in the real part, and of 0{l/a) in the 
imaginary part, it is clear that the double- gradient terms of equation (H) are 
the same as in the free-particle case @. Collecting all the terms discussed 
above, we recover from the general, nonrelativistic Schrodinger equation 

d I (h ^ - \^ 

ih—ip(x,t) = — \ -V - eA(x,t)] i){x,t) + U{x,t)i)ix,t) (14) 
ot 2m \ i J 

with the vector potential of equation (pH]) and the scalar potential 

U(x,t) = Eo + y^Ki(x,n,t) A(x,ty (15) 

^ 2m 

n 

In canonical quantization, the prescription is to identify U{x, t) and A{x, t) 
with the corresponding functions of the classical Hamiltonian. This amounts 
to an extrapolation into the microscopic domain. The corresponding proce- 
dure in the present context is as follows. By Ehrenfest's theorem, eq. (|T^ will 
reproduce the classical equations of motion for the centre of a wave packet 
in the limit ^ — >■ 0. The classical U {x, t) and A{x, t) then coincide with those 
in the Schrodinger equation. Thus, if desired, U{x, t) and A{x, t) can be 
extrapolated as with canonical quantization. 

In concluding the section, it should be noted that the order-of-magnitude 
assumptions for the hopping parameters depend on the further assumption 
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that no dramatic cancellations occur between k{x, n, t) for different n. Of 
course, those cancellations would require some extra reason for a fine-tuning. 
In the absence of a reason, the assumptions describe the most general and, 
hence, the most likely set of parameters consistent with the constraints. 



3 Free electromagnetic field 

This section is to demonstrate that "unitary hopping" can be a useful concept 
also for quantum field theories. We here consider source- free U{1) gauge the- 
ory. Its Hamilton operator in the temporal gauge ^ is an cxD-dimensional 
version of (|T^). A "hopping" scenario requires the configuration space to be 
discrete. Thus local gauge invariance will have to be discretized, too. In case 
of f/(l) this can be done in a way that preserves an exact local gauge group, 
namely Z{N), whose limit N oo reproduces U{1). 

In lattice gauge theory, a gauge field lives on the links between next- 
neighbour lattice sites. A link can be specified by the site s = {nx,ny,nz) 
from which it emanates in a positive direction, and by the corresponding 
/c = 1, 2, 3. In Z{N) gauge theory [j^ the link variables are phase factors of 
the form 

^2.u/N / = o,l,...,Ar- 1 (16) 

They are related to the electromagnetic vector potential A{s,k), integrated 
along the link, by 

exp {2nil/N) = exp (iaeA/h) (17) 
Thus a Z{N) gauge field configuration is determined by the numbers 

l{s, k) = l{nx,ny,nz, k) nj = 0, ±1, ±2, . . . A; = 1,2, 3 (18) 

We shall indicate by omitting the arguments s and k that we mean the 
configuration as a whole. 

The Hamiltonian will be postulated below to be invariant under charge 
conjugation C, and under space inversion V about any point Sq. As it follows 
from the relation ( p!?] ) to the vector potentials (see also 0), C and Vg^ are 
characterized by their action on the link variables, 

Cl{s,k) = -l{s,k) (19) 
VgJ{s,k) = -l{2so - s-k,k) (20) 
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We also postulate invariance under local Z{N) gauge transformations. These 
are characterized by a number g{s) = 0, 1, . . . , A^ — 1 on each lattice site. The 
link field configuration transforms according to 

l'{s,k) = l{s,k)+g{s + k)-g{s) 

The elementary gauge-invariant construct on a time slice is the plaquette 
variable 

p{s, i, k) = l{s, i) + l{s + I, k) - l{s + k, i) - l{s, k) (21) 
Gauge-invariant, too, is any shift of a link variable; in particular, 

l{s, k) l{s, k)±l if and only if l\s, k) l'{s, k)±l 

The gauge field is quantized by assigning a probability amplitude ipil^t) to 
each link-field configuration /. For this "wavefunction" the general form of a 
unitary-hopping equation is 

ih^i){i,t) = Y,<h^i)i'{i + ^ht) (22) 

Al 

Gauge invariance of the process requires, in the notation of (pi]), 

k{1,AI) = k{p, Al) 

Locality of link interactions is not as uniquely defined — a fact being utilized 
with the "improved actions" of numerical lattice gauge theories [10]. We 
shall only consider the simplest realization of locality, assuming 

• Link-changing processes are independent on different links. 

• A plaquette can influence a change on its own links, at most. 

These assumptions correspond to a pre-relativistic, purely spatial notion of 
locality — no reference whatsoever is made to the phenomenon of light. By 
the assumption of independence, a change on k links within the same time 
interval dt will come with a factor of (dt)'^ and will contribute to the time 
derivative in equation (^2]) only for k = 1. Thus the sum over all link-changes 
Al reduces to a sum over one-link changes. For further simplification, we only 
consider a change by one unit, corresponding to nearest-neighbour hopping 
in configuration space. Thus (|2^ ) takes the form 



zh^M,t) = Y.Y.^±ip-^^'')^i^±^s,,t) = Hij{l,t) (23) 

C-'' links ± 
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where 

J 1 on link s, i 
I elsewhere 
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We intend to Taylor-expand the wavefunction. Instead of the derivative d/dl 
on each link we prefer to use the lattice version of the functional derivative 
6/SA with respect to the vector potential. / and A are related through 
equation ([T7|). Hence, d/dl equals the partial derivative {2nh/ eNa) {d / dA) . 
Now d/dA can be expressed by the functional derivative d/6A essentially by 
introducing factors so that in the characteristic relation dA{s, i)/dA{s' , i') = 
6gg'6ii' the 6gg/ is changed into the lattice delta function a~^6gg'. Thus, 

d 2TTha^ 5 



dl{s,i) eN 6A{s,i) 
Expanding the wavefunction up to order a'^ we have 

,p(,±u,.,t) = ^il,t)± — ^^+-^^ (24) 

The first-derivative term is immediately discarded if we postulate that the 
Hamiltonian be invariant under space inversion V (cf. (pOD). This is because 
the plaquette variables in the hopping amplitudes n'^{p\s,i) are invariant 
under V whereas I and hence d/dA changes sign. 

It remains to discuss the multiplicative terms of (0). To expand the 
hopping amplitudes in a power series in a, we note |^ that the magnetic 
flux density Bi = ^eijkFji^ is related to the plaquette variable by 

exp {io?eFjk{s) /fij = exp (27rzp(s, j, k)/N) 

Thus, at a given flux density of C(l), the plaquette phase factor deviates 
from 1 only in (9(a^), while the plaquette variable p is of 0{a'^N). To be on 
the safe side, we therefore expand the hopping amplitude as a function of 
a^Fij instead of p. Furthermore, we invoke our locality postulates to restrict 
plaquettes with an influence on link {s,i) to the four cases p{s,i,j) and 
p{s — j,i,j) with j 7^ i. Thus, expanding n±{p;s,i) to 0{a'^) we obtain 

^^±\s,t) + ^J2{^±i^^hJ)F^M + n^l\s,t,j)F,,{s-j)) + (25) 
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where in the last term we have discarded any shift of s by j or j' as this 
would lead to an 0{a^) contribution. 

The terms of expression (^) must vanish if the Hamiltonian is to be 
invariant under charge conjugation. This is because C (cf. (|1^)) reverses the 
values of both links and plaquettes, hence reverses the sign of the term in 
([25|), while all remaining terms of ( p5D and also of (p^ are C- invariant. 

By translation invariance of the hopping process, all k's must be indepen- 
dent of the site vector s. By invariance under reflections about a coordinate 
plane, K^2\^}h3}j') the term of (^) must be proportional to Sjj^. 
Hence, by cubic rotational invariance, it must be independent of i. For the 
same reason, K,±\i) as the relevant coefficient of 6'^tlj{l,t)/6A{s,i)'^ must be 
independent of i. 

Inserting in (p3|) the remaining terms of (0) and p5|) we identify the 
Hamiltonian of free QED as 
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where v = J2si{i^+^ + i^-^) is the vacuum energy and where 

In the limit a — >^ we put x = as and d?x = to obtain the familiar form 
H = v + - f E\x) d^x + — j B\x) d^x (27) 

2 J 2yUQ J 



where 



EAs) = — 



eo 5A{s,i) 

4 Conclusions 

We have derived the Schrodinger equation for a nonrelativistic scalar particle 
and for the free electromagnetic field, starting out from the superposition 
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principle for state vectors, using the statistical interpretation, and exploiting 
spatial symmetries to a large extent. The ambition was to avoid any use of 
the distinctly non-quantal concept of trajectories, even in the path-integral 
sense. 

In the case of a free particle, which has all the exploitable symmetries, the 
approach taken here should be compared with the general, group-theoretical 



approach to quantum mechanics as exposed, for example, in |[TT[]. The main 
difference is that we found it unnecessary to consider any classical space- 
time symmetries (Galilei or Lorentz transformations). Rather, the structure 
of the dynamics follows from spatial symmetries together with the absence 
of motional information from states such as \x,t). That absence induces 
symmetries of the time evolution which, however, can be realized only by 
way of a superposition. 

As we have seen, Taylor expansions led to 2nd order derivatives and, in 
the case of QED, to the magnetic energy in the Hamiltonian. The sign of 
the Taylor coefficients, though, must be determined by extra arguments. For 
the mass parameter m in equation (0), it is a matter of convention whether 
kinetic energies are always taken as positive or always negative, so both signs 
of m would seem to make physical sense. A similar remark applies to the 
case of free QED, except for the relative sign of the parameters eo and fio in 
(^^. Here an additional assumption is required, such as the existence of a 
ground state, to recover the positive phenomenological sign. 

For the definition of the mass in ([^) it was essential that a free particle 
find identical hopping conditions on every site of the lattice. But this is also 
what characterizes the lattice as a cartesian coordinate system. In case of 
QED, a cartesian structure is comprised in the local Z{N) gauge invariance. 
Thus the unitary-hopping scenario may explain why cartesian coordinates 
play such a preferred role in a wide range of quantum systems 

Within the "physical" subspace of locally gauge-invariant states, the 
Hamiltonian dynamics of the electromagnetic field as described by (|27| ) is 
automatically Lorentz invariant. This is quite remarkable since we derived 
the dynamics from quantum-mechanical principles in which the roles of space 
and time are initially very different. A similar observation was made by 
Bialynicki-Birula |^ with respect to the Weyl equation. 
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